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Abstract 

Kuzmichev recently showed that use of phase velocity rather than group velocity for 
Cherenkov light signals and pulses from calibration lasers in high-energy neutrino 
telescopes leads to errors in track reconstruction and distance measurement. We 
amplify on his remarks and show that errors for four cases of interest to AMANDA, 
IceCube, and RICE (radio Cherenkov detector) are negligibly small. 



1 Introduction 



Kuzmichev [1] pointed out that participants in high-energy neutrino projects 
mistakenly use tabulated values of refractive index, n, which corresponds to 
phase velocity, instead of using values of group refractive index, n g , which 
applies to signals transmitted at group velocity, v g , where n g {uS) = c/v g (uj). 
The choice of n rather than n g or of v rather than v g leads to errors in track 
reconstruction and in distance measurements. 

The subject of which velocity to use for Cherenkov light has an interesting 
history. Tamm [2] was the first to point out (correctly) that for a particle 
moving with speed (3c the wavefront of a Cherenkov "bow wave" is a cone 
(the "Mach cone") with semi-angle 

cot a = (/?V - l) 1 ' 2 + nuj(3 2< ^/{{3 2 n 2 - I) 1 ' 2 (1) 

da; 



whereas the Cherenkov radiation is emitted at the angle 6 given by 

cosfl = -^ (2) 

pn 
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The presence of the second term in eq.(l) results in 9 not being the com- 
plement of a and in the direction of light emission not being exactly normal 
to the wavefront. Unfortunately, most textbooks draw ray vectors normal to 
the bow wave. Surprisingly, even Sommerfeld [3] wrote, "The ray direction is 
perpendicular to the Mach cone." 

Motz and Schiff [4] corrected the Sommerfeld statement, showing clearly that 
"the radiation propagates in directions that make the angle 9 with the path 
of the [moving] charge, where [cos^ = l/(3n], so that the cone may be thought 
of as 'side-slipping' as it moves along with the charge." They derived eq.(l) 
and showed that only if dn/du = would the ray direction be normal to the 
cone. 

We now discuss several ways in which data for AMANDA and other neutrino 
telescopes need to be corrected, and show that the corrections are mostly very 
small. 



2 Direction and velocity of Cherenkov light from a muon in ice or 
water 

The direction of propagation of the Cherenkov light does not need correction, 
but the velocity does. Using an extensive table of data [5] and interpolating 
linearly in log A for phase refractive index for ice, we generated a polynomial 
(with wavelength A in micrometers) 

n(X) = 1.55749 - 1.57988A + 3.99993A 2 - 4.68271A 3 + 2.09354A 4 (3) 

in the region 0.3 to 0.6 /im for which Cherenkov light contributes to the 
signal of a muon. For group refractive index, we generated a polynomial for 
the fractional increase 

= 0.227106 - 0.954648A + 1.42568A 2 - 0.711832A 3 (4) 



where dn/dX is obtained by differentiating (3). The required correction for 
distances traveled by Cherenkov light from points on a muon trajectory to 
phototubes decreases monotonically from -5.0% at the shortest wavelength, 
0.3 /im, to -1.4% at the longest wavelength, 0.6 /im, if n(A) had been used in 
previous analyses. However, the AMANDA collaboration simply used a fixed 
value, n = 1.32, for all wavelengths. The correction to this value would then 
range from -6.1% at 0.3 /im to -0.5% at 0.6 /im. In all cases the true distance 
is less than the value previously used. Taking into account the Cherenkov 
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spectrum, the transparency of the glass pressure vessel, and the phototube re- 
sponse, the most effective wavelength is around 0.38 //m, with a corresponding 
correction of -3% for distances. Phototube jitter is ~10 ns. The strongest cut 
for AMANDA track reconstruction is that N dir > 4, where iV d i r is the number 
of "direct hits," defined as Cherenkov photon arrival times between -15 and 
+50 ns of those expected for no scattering. A 3% error in distance corresponds 
to a 2 ns error in timing, which is much smaller than phototube jitter. 



3 Use of pulsed lasers to measure spacing between strings 

Since distance is inferred from arrival time of the leading edge of a pulse dis- 
tribution, distance is proportional to pulse velocity and inversely proportional 
to refractive index. All horizontal distances within the AMANDA detector, 
i.e. the interstring separations, must be reduced by the factor n/n g . A YAG 
laser with wavelength 0.532 /zm was used for distance measurements. Distances 
must therefore be reduced by the factor 1.32/ra ff (0.532) = 0.99. 



4 Use of pulsed light sources to infer wavelength-dependent ab- 
sorption length, A a , and wavelength-dependent effective scatter- 
ing length, A e 

For distances d 3> A e ~ 25 m, the expression for three-dimensional random 
walk with absorption gives an excellent approximation to the propagation of 
light pulses in ice [6]: 

V ; \Airv g X e tJ 1 \ Av g X e t \ a J y ; 

Ignoring the pre-exponential factor, we see that the estimate of A a scales as 
n g and the estimate of A e scales as d 2 /n g . If distance had been measured with 
a meter stick or other method that did not involve light pulses, the estimate 
of A e would scale only as l/n g . However, we measured distances with light 
pulses, as a consequence of which the net scaling should be as n g . As a specific 
example, for two strings > 100 m apart, eq.(5) is a good approximation, and 
we find that both A a and A e need to be decreased by 1%, for studies carried 
out with a YAG laser. As an independent check, we generated a function 
u(d,t) for values of A a = 20 to 100 m, A e = 15 to 45 m, and distances 100 to 
300 m. We first used the phase refractive index to generate u(d,t); then we 
increased n by fractional amounts ranging from 0.5% to 5% and evaluated the 
new values of A a and A e that gave the best fit. For separations inferred from 
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the pulse method, we found that SX/X increases as (5n/n) K in both cases, 
with re = 1.03 ± 0.02 for A a , and re — 1.01 ± 0.02 for A e . For separations 
assumed to have been measured with a meter stick, the best-fit powers were 
re = 1.03 ± 0.02 for A a and re = -1.01 ± 0.02 for A e . Thus, both methods lead 
to the same conclusion regarding scaling of absorption and scattering lengths. 



5 Use of pulses for timing calibration 

In AMANDA, the individual time offsets (to) between the times recorded by 
the surface electronics and the times when photons hit the phototubes are 
measured with YAG laser pulses (0.532 pm) sent down optical fibers to dif- 
fusive balls near or inside the phototubes. For most phototubes a fiberball 
located within 1 m is used, but for 15% of the phototubes the pulses have 
to travel through at least 10 m of ice. The time for this propagation through 
ice has to be corrected by n g /n = 1.01, which for distances of 10 to 30 m 
corresponds to 0.44 to 1.3 ns. This correction is smaller than the overall mea- 
surement error of individual timing offsets and is considerably smaller than 
the uncertainty in the recorded times of light pulses. 



6 Dependence of group refractive index on density and tempera- 
ture 

AMANDA-B is instrumented at depths from 1.5 to 1.9 km, and IceCube will 
be instrumented at depths from 1.4 to 2.4 km. The change in refractive index 
with density is given by 

dn = (l-A )(n 2 + 2)(n 2 -l) 
dp Qnp 

where Ao is a small correction due to strain-induced polarizability [7]. Over a 
depth from 1.4 to 2.4 km, the density of ice actually decreases slightly, due to 
the increase in temperature with depth. From data of Gow [8], we estimate 
5p/p 0.004 over a 1 km change of depth from 1.4 to 2.4 km, which leads 
to 5n/n ~ 0.002. Over a depth from 2.0 to 2.3 km in seawater (applicable to 
ANTARES), 5n/n ps 0.0004 due to the density increase. 

For IceCube, the ambient ice temperature increases from -41°C to -19°C for a 
change of depth from 1.4 to 2.4 km. The measured change in refractive index 
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of ice with temperature is given [8] by 

^ = -2.7 x l(T 5 n 
dT 



(7) 



which leads to a fractional decrease 8n/n = — 6 x for the 1 km increase 
in depth. In Mediterranean seawater the temperature decreases about 1° per 
km at depths relevant to ANTARES; the resulting change of refractive index 
due to a change of temperature is negligibly small. 



7 Conclusions 

For 0.532 /zm pulses in ice, no error exceeds ~1% when a single value of 1.32 
for refractive index was previously used. If a nitrogen laser (0.337 //m) were 
to be used, the error would grow to 5% relative to the value 1.32. 

For radio Cherenkov radiation in ice (RICE), dn/dX = for wavelengths from 
a few cm to ~10 m, which spans the entire region of interest. Thus, values of 
the phase refractive index can be used without introducing error in analysis 
of RICE data. 

For ANTARES and other neutrino telescopes in seawater, 5n/n w 0.0012 over 
a depth interval of 1 km, due to the dependence of n on pressure; the change 
due to dependence of n on temperature is negligible. 

For AMANDA and IceCube, the change in refractive index due to the increase 
of ice density leads to an error of only ~0.2% per km. The error due to a 
temperature increase with depth would be even smaller, ~0.06%. 
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